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Two new methods for investigation of two-dimensional quantum systems, 
whose Hamiltonians are not amenable to separation of variables, are pro- 
posed. l)The first one - SUSY— separation of variables - is based on the 
intertwining relations of Higher order SUSY Quantum Mechanics (HSUSY 
QM) with supercharges allowing for separation of variables. 2)The second 
one is a generalization of shape invariance. While in one dimension shape 
invariance allows to solve algebraically a class of (exactly solvable) quantum 
problems, its generalization to higher dimensions has not been yet explored. 
Q I Here we provide a formal framework in HSUSY QM for two-dimensional 



quantum mechanical systems for which shape invariance holds. Given the 
knowledge of one eigenvalue and eigenfunction, shape invariance allows to 



' construct a chain of new eigenfunctions and eigenvalues. These methods are 

applied to a two-dimensional quantum system, and partial explicit solvability 
is achieved in the sense that only part of the spectrum is found analytically 
and a limited set of eigenfunctions is constructed explicitly. 
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1. Introduction 



In one-dimensional Quantum Mechanics the importance of exactly solvable and quasi 
exactly solvable (QES) models has been stressed repeatedly. The approach of Supersym- 
metric Quantum Mechanics (SUSY QM) and in particular shape invariance m has been 
fully exploited for construction and investigation of such models by generating a partnership 
between pairs of dynamical systems which allows to establish the solvability of one in terms 
of the other by means of intertwining relations with supercharges of first order in derivatives. 

With this knowledge one can construct a variety of multidimensional solvable quantum 
systems by suitable combination of solvable one-dimensional dynamical systems (separation 
of variables and its generalizations [0). From now on we will focus the main attention 
to two-dimensional quantum systems without assuming that such separation of variables is 
possible. In two-dimensional Quantum Mechanics solvable (or partially solvable) dynamical 
systems, for which the entire spectrum (or part of it) and the associated wave functions are 
known, play a role similar to solvable (or QES) models in one dimension. 

Within the search for a larger class of problems, which can be solved by supersymmet- 
rical methods, extensions of SUSY QM have been elaborated with different realizations of 
the intertwining operators (supercharges). In particular, one-dimensional supercharges were 
constructed in terms of higher-derivative operators 0, (HSUSY QM), the associated 
superalgebra (HSUSY) containing a higher order polynomial of the Hamiltonian. This gen- 
eralization has been revisited recently in 0, and referred to as A^-fold SUSY. Its 
most simplified version {A-tjpe A^-fold SUSY) corresponds to a solvable ansatz for the 
so called 0] reducible supercharges and for more general factorized supercharges. However 
also a class of higher order supertransformations was found 0, which can not be represented 
as a succession of two standard first order supertransformations (the so called, irreducible 
supercharges). Accordingly, the second order intertwining operators were exhaustively clas- 
sified as reducible and irreducible. We warn that our definition of reducibility P] does not 
coincide with factorization of supercharge. 
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An important step for investigation of one-dimensional HSUSY models is the construction 
of zero modes of (higher order) supercharges which can be instrumental for quasi exactly 
solvability and for the study of spectral properties of higher order shape invariant 
systems. Indeed, it is well known that in one-dimensional standard SUSY QM a very elegant 
method of solution of the spectral problem exists for potentials which preserve their shape 
in the SUSY partnership It was shown that for a wide variety of such potentials 

the entire spectrum and the eigenfunctions can be obtained algebraically providing a fresh 
reformulation of the old Factorization Method ^ for the Schrodinger equation. Extensions 
of shape invariance to third order not fully reducible HSUSY QM models were proposed in 
i- 

In the A^— dimensional SUSY QM [|l^] starting from a scalar Schrodinger Hamiltonian 
a chain of matrix Hamiltonians of different dimensionality was constructed ending with 
another scalar Hamiltonian. Each pair of neighbouring Hamiltonians is intertwined by a 
first order supercharges, each Hamiltonian has a partial isospectrality with both neighbours. 
In the two-dimensional case this chain simplifies to two scalar Hamiltonians and one matrix 
2x2 Hamiltonian. The spectra of the two scalar Hamiltonians build up the spectrum of the 
matrix Hamiltonian, but in general the spectra of the scalar Hamiltonians are not related. 

However, cases were found ||12[, where the two scalar two-dimensional Hamiltonians 



are intertwined by second order supercharges and are therefore isospectral, up to zero modes 
of the supercharges. This suggests that HSUSY QM may even be more important for the two- 
dimensional case than in one-dimensional case to study spectra and eigenfunctions of two- 
dimensional models not amenable to separation of variables. Due to the complexity of the 
system of nonlinear partial differential equations, arising from the second order intertwinig 
relations, one has to look only for particular solutions. Indeed, two classes of such particular 
systems were found in ||Tl|, [|T^ . 

In Section 2 we outline the second order SUSY QM framework in one and two dimen- 
sions, which allows to formulate in Section 3 the method of SUSY— separation of vari- 
ables for two-dimensional Hamiltonians, which are not amenable to separation of variables. 
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This method can be used when the supercharges allow for separation of variables. In Sec- 
tion 4 a two-dimensional singular Morse type potential (Subsection 4.1) is shown to satisfy 
Sf/SF— separation of variables. For this model the zero modes of the supercharge (Subsec- 
tion 4.2) and eigenfunctions of the Hamiltonian (Subsection 4.5) in the linear space of zero 
modes of the supercharge are constructed. In addition, a new eigenfunction outside this 
space is built explicitly (Subsection 4.6). The role of shape invariance (Subsection 5.1) has 
not been explored yet in the two-dimensional SUSY QM (Subsection 5.2). It is one aim of 
this paper (Subsection 5.3) to find Hamiltonians (without separation of variables), which re- 
alize the two-dimensional generalization of shape invariance. The model of Subsection 4.1 is 
found to satisfy shape invariance. It is shown that, due to the nontrivial space of zero modes 
of supercharges, shape invariance does not allow a fully algebraic solution for the entire 
spectrum, but only (partial explicit solvability) for part of the spectrum and for the corre- 
sponding wave functions. Section 6 contains a brief discussion of the two-dimensional Morse 
type model in another region of parameter-values and of the integrability of this model, by 
elucidating the action of quantum integrals of motion, commuting with the Hamiltonian, on 
the constructed eigenfunctions. 

2. Second order SUSY in one and two dimensions 

Let us provide notations and the main formulas for standard SUSY QM |[^: 

H = Q+Q- = -d^ + Vix); Mn{x) = ^„^„(x); (1) 

H = Q-Q+ = -d^ + Vix); H^^ix) = ^„^„(x); (2) 

HQ+ = Q+H- Q-H = HQ-; (3) 

Q+ = -d + W{xy, Q- = {Q+)^ = d + W{xy, (4) 

^n(a;) = Q-^nix); ^nix) = g+^„(x), (5) 

where d = d/dx. 
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The spectral equivalence (up to zero modes of Q^) of if, H can be expressed via the 
superalgebra: 





(Q-)^ = 0; [if,Q^] = 0. (6) 



The second order generalization of the superalgebra incorporates |]^ the most general 
intertwining operators of second order in derivatives: 

Q+ = d^ -2f{x)d + h{x)- Q- = {Q+y. (7) 

Intertwining relations with supercharges (0) lead to the expressions of the potentials 
V, V and the supercharges in terms of the only real function f{x) : 



V{x) = -2 fix) + fixY + 



2 , fix) ff{x)Y d 



2 fix) \2fix)J Afix 



^2 , fix) ff'ix)^^ d 



)2 - 

where a and d are arbitrary real constants. The case c? < was called 0] reducible, since 
it can be interpreted in terms of two successive first order supertransformations with real 
superpotentials and an hermitian intermediate Hamiltonian (up to a constant real energy 
shift). The alternative case > was called irreducible. 

A two-dimensional generalization of standard SUSY QM was proposed in , where the 
Superhamiltonian includes a chain of one matrix and two scalar Hamiltonians. Each scalar 
Hamiltonian is separately intertwined with the matrix Hamiltonian, but the spectra of the 
two scalar Hamiltonians are not related. 



If one wants to relate directly |Tl|, [jT2[ the spectra of two scalar two-dimensional 
Schrodinger operators via intertwining relations analogous to (Q): 

Hix)Q+ = Q^Hix), Q-Hix) = Hix)Q-, (8) 
H = -A + Vix); H = -A + Vix) A = dl + d^, di = d/dxi, 



it is expedient to make use of a two-dimensional generalization of (|^: 

= {Q~y = 9rk{mdk + C^{x)d; + 5(f), (9) 

where all coefficient functions are real. 

It means that, up to zero modes of Q^, the spectra of if, H coincide and their eigenfunc- 
tions are mutually connected: 

A classification of the dynamical systems requires first to determine the "metrics" in 
supercharges (^. The Eq.® leads to a set of four equations for the metrics gik{x), 

diQu + d^ga + = 0; {ikl} = {111; 112; 122; 222}, (10) 

which can be solved independently from the other ones: 

Qii = ax\ + aiX2 + hi] g22 = ax\ + 02X1 + 62; 9x2 = -^(2aa;iX2 + aiXi + 02X2) + 63, 

with a, Oj, hi - constants. 

For the particular case of the unit metrics Qik = Sik the corresponding quantum systems 
allow |jTT| for the so called R-separation of variables in parabolic, elliptic or polar coordinate 



systems. 

After having solved (PUD, the intertwining relations (^ are equivalent to the following sys- 
tem of differential equations for potentials V{x), V{x) and coefficient functions Ck{x), B{x): 

d,Ck + dkQ + Ag,k -{V- V)g,k = 0; 

AQ + 2d,B + 2g,kdkV - {V - V)C, = 0; (11) 
AB + QikdAV + CAV -{V- V)B = 0. 

Because of the complexity of the equations above, only particular solutions have been 



found in |11|, |T2| by making suitable anzatzes. In particular, for the supercharges with 



Lorentz metrics {gik = diag{l, —1)) : 

Q+ = (dl - dl) + CA + B = 49+9_ + + + B, (12) 



a solution of (|TT]) can be reduced |TT|, to a solution of the system: 

d^C^F) = -d^C+F); (13) 
dlF = dlF, (14) 

where x± = Xi ± X2 d± = d/dx± and C± depend only on x±, respectively: 
The function F, solution of (p!^, is represented as a sum 



F = Fi(x+ + a;_) + F2{x+ - x^). 

The potentials V{x), V{x) and the function B{x) are expressed in terms of Fi(2xi), ^2(2x2) 
and C±(x±), solutions of system (pTSD, (|14|): 

^ = i(C7; + Ci) + i(C^ + C2) + i(F2(x+-x_)-Fi(x+ + x. 

= -i(C; + Ci) + ^(C^ + C!) + i(F2(x+-x_)-Fi(x++x_)), (15) 
5 = ^(C+C_ + Fi(x+ + x_) + F2(x+-x_)y (16) 



3. /S^J/Sy— separation of variables: a construction of 
new two-dimensional partially solvable models 

We want to study two-dimensional Hamiltonians without separation of variables within 
the SUSY approach. Our goal is to find a class of Hamiltonians for which part of the spectrum 
and of the eigenf unctions can be found (partially solvable systems). Our approach is based on 
the solution of the intertwining relations (P) for partner potentials and supercharges followed 
by the investigation of zero modes of the supercharges (0). Because, in principle, there is no 
connection between separation of variables in the supercharges and the Hamiltonians, one 
can look for the opportunity of partially solving for the spectrum of a Hamiltonian (which 
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does not allow the separation of variables) via normalizable zero modes of supercharges which 
do allow such separation. 

More specifically we will start from the investigation of zero modes of supercharges which 
do allow for separation of variables and solve the intertwining relations (H) obtaining Hamil- 
tonians which are not amenable to separation of variables. The algebraic method (for partial 
solution) will be presented in this Section in its general form. In the next Section the method 
will be applied to a specific type of two-dimensional models. 

Let us suppose that there are (A^ + 1) normalizable zero modes fln{x), n = 0,1, N of 
the supercharge Q~^: 

Q+n{x) = 0, (17) 

where Q{x) is a column vector with components fi„(x). Acting by intertwining (||) onto Q{x) 
it is easy to realize that the space of zero modes is closed under the action of H : 

Hn{x) = Cn{x), (18) 

where C = 11 

^ik I I is cL C number x— independent real matrix. If the matrix C can be diago- 
nalizedQ by a real similarity transformation: 

BC{B)-^ = A = diag{\Q, Ai, A^v), (19) 

the problem reduces to a standard algebraic task within the zero modes space: 

H{M{x)) = A{M{x)). (20) 

For attacking the problem it is expedient to eliminate the first order derivative terms in 
the supercharge (|12D by a suitable similarity transformation: 

g+ = exp{-xmQ^ exp{+x{x))=dl-dl + ^{F,{2x,)+F2{2x2)), (21) 

X{x) = -^{Jc4x+)dx+ + Jc4x-)dx^y (22) 

''if C can not be diagonalized, the algebraic method can be apphed if one can solve the less restrictive 
matrix equation BC = AB. 
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We notice that exhibits separation of variables: this is what we mean by 
Sf/SF— separation of variables, even if the components of the Superhamiltonian do not 
admit such separation. Zero modes of can be found as a linear superposition of products 
of one dimensional wave functions rjn{xi), and p„,(x2), satisfying Schrodinger equations: 

(-c^i - ^-Fi(2a;i))?7„(a;i) = e„r7n(a;i), 

{-dl + ^F2{2x2))pn{x2) = e„p„(x2), (23) 

with e„ - the separation constants. 

In analogy to (|2l|) , one can define operators 

/, = exp(-x(f))// exp(+x(x)) = -9?-92'+Ci(f)9i-C2(x)92-^Fi(2xi) + ^F2(2x2), (24) 
h and eigenfunctions of g"*": 

ujn{x) = exp(-x(x)) ■ ^ln{x), (25) 

keeping however in mind that the normalizability and orthogonality is not preserved au- 
tomatically due to non-unitarity of the similarity transformation. The operators h and h 
should not be interpreted literally as Hamiltonians since they are non-hermitian (but have 
real spectrum) 0, and are not factorized in q"^, like in (|l]), (||). 
Then using (P3| ) one can write: 

hUnix) = [2en + Ci{x)di - C2ix)d2]cOnix). (26) 

While from the equation above it is not manifest, we know however from (|T8|) that the space 

spanned by functions u!n{x) is closed under the action of h. In the concrete model of next 

Section this will be demonstrated in Subsection 4.4. 

"^Incidentally, we remark that by an additional unitary transformation one can generate analytically a 
class of two-dimensional non-trivial Hamiltonians with complex potentials but with real spectrum. Examples 
of two-dimensional Hamiltonians have been discussed recently in . 
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It is clear that, in contrast to (^) where variables are separated, there is no separation 
of variables for /i, which makes the two-dimensional dynamics not-trivially reducible to one- 
dimensional dynamics. For this reason we refer to this method for partial solvability as to 
Sf/SF— separation of variables. 



4. A new two-dimensional partially solvable model. 
4.1. Formulation of the model. 

In this Section we apply the method of Section 3 to a particular example, which can 
be interpreted as a suitable two-dimensional generalization of Morse potential. Among the 
solutions of the system of equations ( |T3| ) and (|13|) found in [|I2| we focus attention on the 
particular case already presented in |T^: 

Fi{x) = /ci ^a;+a_ exp(Ax) -|- exp(— Aa;)^ + A;2 ^tt+tt^ exp(2Aa;) + exp(— 2Aa;)^ ; 

—^2(0;) = ki exp(Aa;) + exp(— Ax)^ + k2 exp(2Ax) + exp(— 2Aa;)^ ; 

^ _^ a±exp{Xx±) + I3±exp{-Xx±) ^^y^ 



A^a±exp(Ax±) — /3± exp(— Ax±)^ 



^ _ 2a+p+{l + 8A2) + al exp(2Ax+) + Pl exp(-2Aa;+) ^ 
SA^ ^a;+ exp(Ax+) — /3+ exp(— Aa;+)^ 
2a_/5_(l - SA^) + al exp(2Ax_) + pi exp(-2Aa;_) 
8A^ a_ exp(Aa;_) — /5_ exp(— As_' 



ki a+P- exp(2Aa;2) + exp(— 2Aa;2) 1+^2 exp(4Ax2) + 



1 

+a1p'^ exp(— 4Aa;2) j + ^1 {^a+a- exp(2Axi) + P+P- exp(— 2Axi) ) + 
+fc2( a+a-exp(4Axi) + exp(-4Aa;i; 



(28) 



We will consider a specific case of the expressions for C± and V by choosing the parameters 
in (p7|) and ( pH]) so that /5+ = 0, a_ = A = —a/ 2, a > 0. It is evident that Eq.(|13|) admits 
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a more general solution than ( pT] ) because one can introduce an additional multiplicative 
parameter in C± in (pT]). The nonlinear dependence of (pIS]) on C±{x±) leads to a nontrivial 
generalization of the corresponding potential. 
We then obtain: 

C+ = Aaa; C_ = 4aa ■ coth (29) 
= -Fi(2xi) = — A^exp(— 2«xi) — 2exp(— axi)^; (30) 
f2{x2) = ^F2{2x2) = +A(^exp{-2ax2)-2exp{-ax2)y, (31) 
. A(2.-l,.„h-(^)+4aV + 

+ A exp(— 2q;xi) — 2 exp(— axi) + exp(— 2ax2) — 2 exp(— 0x2) ; 
V{x) = a^a(2a + 1) sinh~^(^^^^ + 4a^a^ + 

+ A exp(— 2axi) — 2exp(— axi) + exp(— 2^X2) — 2exp(— ax2) , (32) 

where A is an arbitrary positive constant, and a is a parameter originating from the new 
multiplicative constant mentioned above. We will show in the next Sections that the range 
of variation of this parameter a will characterize the dynamics of the model. It is perhaps 
interesting to remark that the reflection a —>■ —a signals the supertransformation Q'^ ^ 
Q~ and H ^ H . Both potentials V{x) and V{x) are also invariant under the interchange 
Xi ^ X2 ("x_-parity" conservation). As for standard P-symmetry, this invariance leads to 
classification of eigenfunctions according to their "a;_-parity" values. 

One easily recognizes in (^) a sum of two Morse potentials plus a hyperbolic singular 
term which prevents to apply the method of separation of variables for the system (|32D. These 
singular terms in V{x),V{x) can be both attractive, for the case \a\ > i, or one repulsive 
and one attractive, for the case \a\ < |. The parameter a will be further constrained by the 
condition that the strength of the attractive singularity at a;_ ^ should not exceed the 
well known bound — l/(4x^ ). 
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4.2. Zero modes of 

The reason for the choice of the model (^) - (|32D is that the corresponding equations 
for the zero modes (^) of can be solved analytically in this case. In particular for 
the discrete spectrum e„ < the normalizable eigenf unctions are: 

Unix) = exp(-^i^)(ei6)'"i^(-n, 2s„ + 1; ii)F{-n, 2s„ + 1; 6), (33) 

where F{—n,2sn + 1;^ is the standard degenerate (confluent) hypergeometric function, 
reducing to a polynomial for integer n, and 



2\'A 

= exp(-axi); (34) 



a 



A 1 

n - - > 0; 35 

a 2 

a , 1 

1 ^(n + - 

"A" 2 



-A 



(36) 



The number (A^ + 1) of normalizable zero modes (|33| ) is determined by the inequality (|35|). 
As a final remark, let us point out that the zero modes iOn{x) of are entirely based on the 
"Morse-part" of the potentials, and consequently do not depend on the parameter a, which 
refiects the strength of the singular part. 

4.3. Normalizability of zero modes of 

We want to discuss the normalizability of the zero modes fin(^) which are connected 
with (H): 

fi„(f) = exp(x(x))u;„(f). (37) 

While in principle one can investigate the most general conditions for which f2„(2^) is nor- 
malizable but iOn{x) is not, we will restrict ourselves for simplicity to the conditions for 
which the normalizability of ujJ^x) implies that for ^^[x). A constructive discussion of these 
restrictions will be given in the following. 

/^From (|22|) and ( p9l) one can find the analytical expression: 



exp(x(x)) = exp(-a«x+)| sinh( ^^jl = (^-^ ■ 1^:^^^^ ) ; a > 0. (38) 
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For ^1 — > ^2 the term |^2 — ^il^^'^ in ( pSD requires for the normahzabihty of ( pT] ) that a < 1/4. 

At infinity (in the > 0,^2 > quadrant) the (|38D does not change essentially the 
behaviour of ujn{x), which is exp( — ^^^^^ ). 

At the origin (again in ^1,^2) the behaviour of u!n{x) ~ (^1^2)'^", s„ > combines with 
behaviour at the origin of (pH). In polar coordinates 



6 = ^ ■ cos 0; 6 = ^ ■ sin > 
the relevant integrand of ||i7(af)|p at the origin reads: 



2{2a+s„) 

sin(2(/)) ' 



e sin(20) |sin(0-|)|4a 

Combining all these restrictions, normahzabihty of fi„(a;) therefore is ensured by: 

1 

s„ + a > 0; Sn + 2a > 0; s„ > 0; < ^• 

In addition, one has to impose the constraint that the singularity for the Superhamilto- 
nian (i.e. both for H and if) should be repulsive or, if attractive, should be bounded, as 
explained at the end of Subsection 4.1, by 4a(2a ± 1) > —1/4. All these constraints can be 
implemented contextually leading to three parameter families of models: 



1 1 , VA 1 

a G (—00, ^): Sn = n > —2a > 0. (39) 

^ ' 4 4^2^ a 2 ^ ' 

Inequalities ( PU[ ) can be satisfied by appropriate choice of parameters a and A, and/or by 
suitable restriction of the considered number N of zero modes Qn{x)- It is perhaps interesting 
to remark that the reflection symmetry H ^ H, ^ for a —>■ —a is broken by the 
fact that the operator Q~ has no normalizable zero modes when the parameters range in 
(^), while has. An alternative possible range of parameters which allows to implement 
explicitly the reflection symmetry will be discussed in Section 6. 

4.4. Algebraic solution. 

In this Subsection we will prove explicitly that the operator h, when acting in the linear 
space of the zero modes of q~^, leaves this space invariant. Inserting (|33| ) into expression ( p6D 
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one obtains: 



hujn{x) = -2{2aa Sn-e„)u;„(f)+- — r'77 rv'^^P(~o 2s„+2;^i,^2) 

ZSn + i (,^2 — ^ 

(40) 

where the function $(6, c; ^1,^2) is defined by: 

Hb, c- ei, 6) = [F{b, c; ^,)F{b - 1, c - 1; 6) - c; - 1, c - 1; ^i)]- 

Relations between contiguous hypergeometric functions lead to: 

m c; 6,6) = - ^i)[F{b, c; ei)F(6, c + 1; 6) + 

^""^ ^ ei6$(fe+l,c + 2;ei,6). (41) 



(c-l)c c(c+l) 

By making repeated use of (^) for 6 = 1 — n one arrives after laborious manipulations 

to: 

^l-n,c■, 6,6) = (6-6)S^=danfc(66)'-[i^(l-^+fc,c+2fc+l;6)i^(l-^+fc,c+2fc+l;6)], 

(42) 

where 

n + c-1 

a-no = — 7 7T— ; a„fe = for k > n; 

(c- l)c 

(n — 1)! (n + c — l)(n + c)...(?T, + c + n — A; — 2) 

""'^ " ^ ■ [(c - l)(c + 2n - 2A; - 2)] • [c(c + l)...{c + 2n - 2k - 3)]^ ' ^ 

In our case c = 2s„ + 2 and (^) becomes: 

$(1-72, 2s„ + 2; 6, 6) = (6 - 6)S^=d«nfc(66)"-'-' i^(-fc, 2sfc + 1; 6)i^(-fc, + 1; 6)- 

Inserting (^) into (|iD| ) and taking into account that Sk — Sn = n — k (see (^^), an equation 
of the type ([TSl) is obtained: 

4 ^ 

/icc'„(f) = Sf^oCnfct^fc(^) = -2(2aa^s„ - e„)t^„(f) + ^ -^ ^IZoank^^kix), (43) 

showing that the matrix C is triangular. We recall (see (|37|)) that the zero modes f2n(^) are 
related to ujn{x) by the similarity transformation. 
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4.5. Eigenfunctions of H. 

The triangular matrix C with all different and nonzero diagonal elements can be diago- 
nalized by the similarity transformation (|19D and its eigenvalues coincide with the diagonal 
elements c^k- From (^) one can construct a set of eigenfunctions ipn{x) and \E'n(x) of h and 
H, provided matrix B is obtained (see Eqs. ([T9|) , (^)): 



The index giving the numeration of the wave functions ip will be elucidated below. The 
eigenvalues of H which correspond to \E'fc(x) are expressed in terms of the parameters of the 
problem: 

Ek = Ckk = ^N-k = -2{2aa^Sk - et). (45) 

lYiom. (^) we conclude that assumptions made for the diagonal elements of C are not very 
restrictive due to the interplay of the parameters of the problem. 
We start from the formal solution for the matrix elements of B : 

h 



(46) 

N—m,p 



where the (A^ + 1) triangular matrices labelled t^^\ m = 0,1,. ..,A^ are defined via the 
matrix elements of C : 

(m) Cn^k 

'^n,k = • 
CN-m,N-m ~ Ck,k 

We stress that in ( ^61) the expression ^r'^™'-'^ means the /-th power of the matrix t^'^\ The 
repeated index (A^ — m) is frozen in (^) and not summed over. This expression allows to 
write all elements of the m-th row hm,p in terms of the matrix r'-'"^ and the arbitrary value 
of the element hm,N-m on the crossed diagonal. These arbitrary values can be fixed by the 
normalization condition for the wave functions N-n{S) in O)- 

/^From the triangularity of C and r(™) in iM] it follows that non-zero elements hmn are 
obtained only for m + p < N. This means the matrix B vanishes below the crossed diagonal. 
It is now clear that the index of the wave function in (|4^ ) is taken in a way to make \l/fc(x) 
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a linear combination of the first {k + 1) zero modes ^li{x); I = 0,1,..., A;. In particular, 
^o(^) ~ ^o(^)- 

We sketch a constructive algorithm to justify the result (^Bf) for the matrix B. It is not 
difficult to see that indeed the system of equations 

'^k=0^i,kCk,l = Kbi,l 

is solved for i = 0, / = by (|46|) with Aq = cat^tv and arbitrary 6o,Af 7^ due to the 
triangularity of C. One can also check that the solution holds for i = and I = (A^ — 1), (A^ — 
2),...,0 by solving iteratively the corresponding linear equations. For i = 1 one convinces 
oneself that the first equation with our assumptions implies bi ^ = 0, which signals the 
crossed triangularity mentioned above. Following the steps as before / = (A^— 1), (A^— 2), 
one again can solve iteratively with arbitrary element 6i,(Ar_i). Similar manipulations can be 
performed for higher values i = 2,3..., N. 

All these wave functions \E'fc(x) live in the space of zero modes, in the next Subsection 
we will construct also additional eigenfunctions, not linear combinations of zero modes. 

4.6. Additional eigenfunctions of H. 

The eigenfunctions of the previous Subsection may be used for constructing more general 
eigenfunctions of h and of H via a product ansatz: 

(f){x) = Mx) ■ Q{xy, ^x) = ^o{x) -Qix). (47) 

The eigenvalue equation for h leads to: 

Le{x) = 76 (f), (48) 

where 

L = -a'^dl + - 2so - 1)% - 4aa'^^{d^, - %). 

The sum over index i is implicit. More specifically, the eigenvalue equation 

h(j){x) = Coo0(x) + i'o{x)L(c>{x) 
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determines the new eigenvalues 



(49) 



For the spectral problem (^) we can provide only particular solutions by choosing suit- 
able ansatzes for the functions 0(x). It is useful to change variables: 

1 1 

so that 



1 1 



Z2 



L = - — 
2 



Z2 



The action of L on monomial products z^' ■ is: 



T 



(r{Pi,P2) + 4(/^2 - 1 - 4a)/52 + 4a(A - 1) + — 

Z2 Zl Zl 



(z'l'-zh, (50) 



where constants (y{Pi,P2) read: 

a{f3i, (32) = A(A - 1) + (32{(32 - 1) + - 4a(2A + /^s) - 2(2so - 1)(A + /52). 

for different powers {Pi, P2) and 



Using ( |5nD and linear combinations of two terms 
{^1,^2), one can construct only three solutions of 

1) 0«(f 

2) e(')(^ 



^(4a+l). 

^2 ) 



7 



(1) 



4a - 2so + 3 



7 



(2) 



3) e^^\x) = zi 



7 



(3) 



a2(2so-l)(4a + l); (51) 
4a\so - l){2a + 1); (52) 

a2r4a + 2so-iy (53) 



4a + 2so - 1 ' 

The eigenvalues 7 of (|48|) can be easily identified, when the corresponding functions G(x) 
do not spoil the normalizability of \E'o(a;). Within the bounds imposed ( pUj ) only the case 
( |53|) is acceptable provided one also imposes sq > —a + |. 

With a similar procedure one can try to construct further eigenfunctions of H based on 
other \E'fc(x); A; > in turn formed via superpositions (^4]) of zero modes. Of course, this task 
will be more difficult than the one we illustrated above, since the coefficient functions of the 
operator L will contain explicitly the hypergeometric functions from (|33|). In the particular 
case n = 0, described above, the hypergeometric functions reduce to 1. 
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5. Shape invariance for two-dimensional systems. 



In Section 4 we studied the two dimensional spectral problem starting from the linear 
space spanned by zero modes f2„(x) of supercharge and then constructing an additional 
eigenstate (|53D outside this space by choosing a suitable ansatz. It is obvious that algebraic 
methods to extend this construction to larger space of eigenf unctions are highly welcome. In 
this respect the well known (in one dimension) method of shape invariance (Subsection 5.1) 
is very suggestive. While in one dimension shape invariance amounts effectively to exact 
solvability, in two dimensions we will show that one can achieve partial solvability only. The 
reasons of this will be particularly clear from an analysis of simple two-dimensional systems 
with separation of variables (Subsection 5.2). On one side, one knows how to solve the 
problem for this system by standard shape invariance for each degree of freedom. On the 
other side, we also will consider this problem directly from a two-dimensional point of view 
and show that, in general, only partial solvability will be achieved in this last approach. 
Finally (Subsection 5.3), we will extend our method to the two-dimensional systems with 
Sf/^y— separation of variables, already described in Section 4, where only partial solvability 
holds. 

5.1. One-dimensional shape invariance and solvability. 

For reader's convenience we write the basic steps of standard shape invariance in absence 
of spontaneous SUSY breaking ||l|] , . One refers to shape invariance when a one-dimensional 
Superhamiltonian H depends on a parameter a and, in addition, its components H and H 
satisfy: 



where d = d{a) is some new value of parameter, which depends on a, and TZ{a) is a (c- 
number) function of a. The absence of spontaneous breaking of supersymmetry for all values 
of a implies that the lowest eigenvalue E^la) of H{a) vanishes and the corresponding eigen- 
functions '^o{(^) are normalizable zero modes of Q~^{a). It is well known |Q, [||, that the 




(54) 
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intertwining relations 

Q-H{a) = H{a)Q^{a) (55) 

with the standard first order supercharge @ allow in this case to solve the entire spectral 
problem for H{a). 

The crucial steps follows. Start from 



Consider the relation fBl) to obtain: 



(56) 



^(a)^o(a) = 7^(a)^o(a)• (57) 
It is important to remark that \E'o(a) = ^o('^) has no nodes and therefore is the ground state 



wave function of H(a). The combination of ([55| ) and (|57D yields: 



H{a) 



Q (a)^o(a) 



7^(a) 



Q (a)^o(a) 



(58) 



Provided 



g-(a)^o(a) 



is normalizable, we have generated an excited state of H{a), and 
thus TZ{a) is naturally positive. It is clear that these steps can be repeated up to the last 
step, where the resulting wave function \E' will no more be normalizable. There are notorious 
cases (oscillator-like potentials) where the spectrum is not bounded from above. 

It is also clear that the isospectrality of H{a) and H{a) (up to the only zero mode \E'o('^)) 
implies that there is no eigenvalue of H{a) between zero and the ground state energy Eo{a) 
of H. This observation leads to a proof that after suitable iterations one gets the entire 
spectrum of H{a). This method is referred as algebraic solvability (or complete solvability) 
by shape invariance in one-dimensional SUSY QM. 



5.2. Two-dimensional shape invariance for systems with sepa- 
ration of variables: solvability or partial solvability. 

Already the trivial two-dimensional model with separation of variables: 

H{x) = Hi{xi) + H2{x2); x = {xi,X2) 
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shows that there is a considerable difference with respect to the one-dimensional case. The 
crucial reason is that the space of zero modes of supercharges becomes now of higher di- 
mensionality including the products of one-dimensional zero modes of the first Hamiltonian 
times all states of the second Hamiltonian and vice versa. 

In order to realize a nontrivial intertwining relations one can consider factorized super- 
charges of second order written as products of first order supercharges^ 

Q^ = Qi-Qt; Qt = {Td^ + w,{x,)). (59) 

Now suppose that Hi and H2 both are shape invariant: 

iJi(ai) = Hi{ai) + 7^l(al); #2(02) = ^^2(02) + 7^2(a2), 
i.e. H is shape invariant with a vector parameter a = (ai, 02) : 

if (a) = //(a) +7^(a). 

While iterations analogous to the previous Subsection are obviously possible, it is clear 
that one can not argue about the entire solvability of the spectral problem, because in general 
many zero modes of (^) exist. Their number depends on the confining properties of Hi and 
H2. For example, in a case of oscillator-like potentials this number becomes infinite, and 
they are distributed over the whole spectrum. In this case only partial solvability of H 
can be achieved by the choice of (|59D and shape invariance. Of course, one can solve such 
trivial models by separate use of = Qf, which allows to solve the entire spectrum of 
two-dimensional model in terms of the one-dimensional ones. 

■^Note that in such models H allows the separation of variables, but do not. If one wants also to 
allow the separation of variables, then one would have to consider Q"^ — Qf. 
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5.3. Shape invariance and partial solvability for two-dimensional 
systems. 

Let us suppose to have two-dimensional system with a Hamiltonian if, which is related 
to H by (|5^) . For simplicity, we have assumed^ that shape invariance is realized only with 
the parameter a. Two-dimensional SUSY QM, realized via ® - @, does not identify zero 
modes of with the ground state of the Hamiltonian. Thus one has to repeat the steps 
- (|^) of Subsection 5.1 by taking into account Eq^o) ^ 0. In order to make our discussion 
more explicit we will from now on refer explicitly to the model ([29| ) - (|3^) with the parameter 
a being bound to (0), as described in the Section 4. 

First of all we observe that this model is indeed shape invariant with 



a = a — -; IZ^a) = Q;^(4a — 1). 



(60) 



We remark also that the infinite domain given by (|39| ) allows iterations of (|5^). The starting 



point is to write (|58 
H(a 



Q-{a)^oia--) 



Eo{a--)+n{a] 



Q (a)^o(a- 



(61) 



where EQ{a) and \E'o(a), not to be identified with ground state, are given by (|4^) and 
Thus we have constructed the new eigenstate and eigenvalue of H{a), provided Q~{a)'^Q{a — 
|) is normalizable. We notice that the eigenvalue I Eo{a — h) + TZ{a) ) is larger than Eo(a) 



with the bounds of (p^). 

It is interesting to compare this first iteration ([61| ) of shape invariance with the solution 
(^) obtained in the framework of the ansatz described in Subsection 4.6. Their eigenvalues 
coincide precisely : 



E = Eo{a)+j 



(3), 



Eo{a --)+ n{a) = a^[Aa{l - so) + (2so - 1)] + 2ec 



while the eigenfunctions (both vanishing for a;_ = 0) differ by a factor reflecting the 

opposite values of "x_-parity". 

^We recall that, in general, there is no connection between the dimensionality of the Schrodinger equation 
and the dimensionality of the parameter manifold. 
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The next iteration of shape invariance will give: 
Hia 



g-(a)Q-(a-^)^o(a-i: 



Eo(a-l)+7^(a-^)+7^(a; 



Q-{a)Q-{a-]^)^^{a-l] 



(62) 

and the new eigenfunction Q~ {a)Q~ {a — ^)\E'o(a. — 1) can be written explicitly as function 
of X. Provided normalizability is ensured, one can thereby construct a chain by successive 
iterations of (|6T|) and (|62D, since Q~{a) has no normalizable zero modes in (^). The end 
point of such a chain will be given by non-normalizability of the relevant wave function. 

6. Discussions and conclusions. 

We want to point out that Sf/SF— separation of variables can be implemented com- 
pletely independently from shape invariance. To this end we present a model where 
Sf/SF— separation of variables holds but shape invariance does not apply. 



Consider the model ( pof ) - (|3^ ) with a choice for parameter a alternative to 



aG(-ii); So>2(|a| + l). (63) 

It is obvious that shape invariance does not apply since the domain ( |63D is too small. One 
important property of this model is that the reflection symmetry a — > —a is explicitly 
implemented. This allows normalizability of both zero modesQ VL{a) = Q{—a) of Q~ and 
f2(a) of and also of all three wave functions $^*^(x) with Q^^\x) given in (^TJ) - (|53|). 

Although we have achieved only partial solvability, we now present a short discussion 
of the quantum integrals of motion (symmetry operators), which exist 0, |Tl[], [|I^ for 



all HSUSY QM systems. Indeed, the intertwining relations (^ lead to existence of the 
symmetry operators R, R for the Hamiltonians H, H, correspondingly: 

R = Q+Q-, R = Q-Q+, [R,H]=0, [R,H] = 0. (64) 



''We recall that uJn{x) by construction do not depend on the parameter a (see Subsection 4.2). 
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In 1-dim case these quantum integrals of motion i?, R become polynomials of H, H with 
constant coefficients. The distinguishing peculiarity of two-dimensional case is given by [|lT 



nontrivial symmetry operators i?, R which are not reduced to functions of the Hamilto- 



nians if, H, i.e. all two-dimensional systems ([15|) which solve the intertwining relations 
are integrable. 



More specifically, it was shown in in the case of second order HSUSY, that for the 
particular case of the unit metrics gik = Sik the fourth order operators R, R can be written as 
second order differential operators up to a function of H,H (R-separation of variables 0). 
For all other metrics Qik the operators R, R are of fourth order in derivatives. 

By construction, the quantum integral of motion R = Q^Q^ gives zero when acting 
onto the eigenf unctions of Hamiltonian from (|4^) . A direct calculation shows that the 
additional eigenfunctions from (|5lD , (0), (153D are also eigenfunctions of operator R : 

with corresponding eigenvalues: 

ri = -(4a + l)(2so - l)(4a - 2so + 1), 
r2 = -16(2a + l)(so - l)(2a - sq + 1), 
ra = -(4a - l)(2so - l)(4a + 2so - 1). 

Thus the eigenfunctions (0) and (|5T| ) - (|53|) are part of the system of common eigenfunc- 
tions of two hermitian mutually commuting operators, the Hamiltonian H and the quantum 
integral of motion R. 

In conclusion, we have formulated two new approaches to partial solvability of two- 
dimensional quantum systems: 

1) The S't/S'F— separation of variables method (Section 3) is to be considered as a par- 
ticular branch of the method of separation of variables, though not for the Hamiltonian, but 
for the supercharge. It is in contrast to the more trivial case, where H allows separation, 
but is factorized and does not allow for the separation (Subsection 5.2). 
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2) The method of shape invariance, well known for the one-dimensional SUSY QM, has 
been reformulated (Subsection 5.3) for the two-dimensional systems, which depend para- 
metrically on a, in the same way as for one dimensional shape invariant systems. This 
construction is based on the knowledge of the eigenfunctions \l'o(ct) and of the eigenvalues 
Eq^g) (for a in the domain (0)). We have illustrated this method applying (Subsection 5.3) 
it to a "singular" two-dimensional Morse system (Subsection 4.1) with Higher order SUSY 
QM in presence of a variety of zero modes of Q'^. 
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